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Abstract. We analyze data on pulsation amplitudes and phases for two (5 Cephei stars, 8 Cet and v Eri. Stromgren 
photometry and radial velocity measurements are used simultaneously to obtain constraints on mean parameters 
of the stars and identification of the excited modes. The inference about the radial mode order and mean star 
parameters is based on comparison of certain complex parameter, /, determined from data, with its counterpart 
derived from linear nonadiabatic modelling of stellar oscillations. The theoretical / values are very sensitive 
to the adopted opacity data. In our modelling we rely on the data from OPAL and OP projects. Significant 
differences were found. New seismic models of v Eri were constructed with both the OPAL and OP opacities. 
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1. Introduction 

In our earlier papers (Daszyhska-Daszkiewicz, Dziem- 
bowski & Pamyatnykh 2003, Daszyhska-Daszkiewicz et al. 
2004) we introduced a new asteroseismic probe of stellar 
structure. The probe is the complex parameter, /, describ- 
ing the ratio of the relative luminosity variation to the 
relative radial displacement of the surface. Empirical / 
values may be derived from pulsation data upon adopting 
an appropriate stellar atmosphere model. The correspond- 
ing theoretical values come from linear nonadiabatic cal- 
culations of stellar pulsation. Comparison of the empirical 
and theoretical / yields a probe of subphotospheric layer. 
Specifically, this is the layer where the mode period is 
comparable to the thermal time scale. It means that the 
layer is not too shallow but still located mostly above the 
frequency-forming interior. Thus, it is poorly probed by 
the frequency data. 

Our new diagnostic tool was first applied to a sample 
of 8 Scuti variables (see references above). In that work 
we were mostly interested in probing the efficiency of sub- 
photospheric convection, whose treatment is a notoriously 
difficult aspect of stellar interior physics. Calculated / val- 
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ues proved very sensitive to the parametrization of the 
mixing-length formalism. 

In B-stars, treatment of convective transport is not a 
significant factor in the modelling uncertainty but it is 
still of interest to ask how well the empirical and theoret- 
ical / values agree. The driving effect in B-type pulsators, 
such as j3 Cephei and SPB stars, arises in the layer of 
the iron opacity bump. This layer has also the dominant 
influence on the / parameter. What is essential for the 
pulsational instability is the location of the layer and the 
iron abundance. Location depends primarily on the effec- 
tive temperature. If the overall metallicity parameter, Z , 
describes the iron abundance in the opacity bump layer, 
then the two crucial parameters, both for the instability 
and the / value, are T e ff and Z. In principle, these two pa- 
rameters can be derived from spectroscopy. However, there 
are considerable uncertainties in the derived values which 
affect accuracy of asteroseismic sounding. We have seen 
this in the case of the v Eri star (Pamyatnykh, Handler 
& Dziembowski 2004, Ausseloos et al. 2004), where the 
constraints on the convective overshooting from the core 
were found to be strongly dependent on T c g. In addi- 
tion, we cannot be sure that the iron abundance in the 
driving layer is the same as in the photosphere. In fact, 
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Pamyatnykh et al. (2004), facing problems with excitation 
of certain modes in v Eri, suggested that iron may be ac- 
cumulated around the opacity bump. Such an effect would 
surely leave an imprint on the value of /. Perhaps such an 
accumulation could also explain the unexpected discov- 
ery of a large number of B-type pulsators in Magellanic 
Clouds (Kolaczkowski et al. 2004). 

It is always important to use all available observables 
to check the accuracy of stellar modelling. A comparison 
of empirical and theoretical / values provides a stringent 
test of models of the atmosphere, envelope and pulsation. 
In addition, the comparison offers an opportunity to test 
microscopic input physics, especially the opacity data. 

Determination of / is possible only in conjunction with 
identification of the I degree for the excited mode. The val- 
ues of I are important for association of measured frequen- 
cies with normal modes of stellar oscillations. There are 
several methods of I identifications based on photometry 
or on spectroscopy, but they required a priori knowledge of 
/ and, frequently, the answer is not unique. In this paper 
we present the first application of our method of simul- 
taneous determination of I and / from combined multi- 
colour photometry and radial velocity data to (3 Cephei 
stars. 

In the next section we recall our method for simulta- 
neous determination of the / parameter and the t degree 
of an observed mode from multicolour photometry and ra- 
dial velocity data. Section 3 is devoted to the analysis and 
interpretation of observational data on 5 Ceti star, which 
is a monoperiodic (3 Cep variable. The analysis of a mul- 
tiperiodic (3 Cep star v Eridani is a subject of Section 4. 
A summary and discussion are given in Section 5. 

2. The method 

We repeat here the basic formulae underlying our method. 
The local radial displacement of the surface element is 
written in the standard form, 

Sr(R,6,(fi) = RReieY^e-' 1 ^}, 

where s is a small complex parameter fixing mode ampli- 
tude and phase, whereas 



YT{9,*) = (-1)*^ {U + W J,™')' (cos 0)e^. 

With the adopted, non-standard, normalization of spheri- 
cal harmonics, \e\ is the r.m.s. value of Sr/R over the star 
surface. The corresponding changes of the bolometric flux, 
JtjoI, an d the local gravity, g, are given by 



and 



- Re{efY e m e- iut }, 
uj 2 R 3 \ Sr 



g \ GM J R 



The complex parameter, /, which is the central quantity 
in this work, has been already defined in the Introduction. 



Since both e and / may be regarded as constant in 
the atmosphere, we can use the static plane-parallel ap- 
proximation. Then, the complex amplitude of the relative 
monochromatic flux variation can be expressed as follows 
(see e.g. Daszyhska-Daszkiewicz et al. 2002): 



where 



A x (i) = V}(if) + £*s, 



e = eY e m (i,0), 



(1) 



V} = b: 



(2 + £)(!- 1) 



l dlog(T x \bj 
4 dlogT eff 

u 2 R 3 n 
+ 2 



GM 



d\og(T x \b^ 



d\ogg 



and 



h x (n)nPt(n)d^. 



The A index identifies the passband. The partial deriva- 
tives of J-\\bg | may be calculated numerically from tabu- 
lar data. In this paper we use data from Kurucz's models 
(2004). For the limb darkening coefficients we take Claret's 
fits (2000). To convert the amplitudes to magnitudes, the 
right hand side of Eq.(l) must be multiplied by the factor 
(-1.086). 

A set of of observational equations for complex un- 
knowns (if) and i is obtained from Eqs.(l) written for a 
number of passbands, A . On the left-hand side we have 
measured amplitudes, A x , expressed in the complex form. 
If we have data on spectral line variations, the set of equa- 
tions (1) may be supplemented with an expression relating 
i to complex amplitudes of the first moments, A4± , 



M\ = \ujR Up + 



GMv* 
R 3 u 2 



e. 



(2) 



where 



h\(n)n 2 Pe(n)dn. 



Jo 

Another useful quantity from spectroscopy is the com- 
plex amplitude of equivalent width variations. This yields 
an additional constraint on / and £. We plan to implement 
this in future work. 

The set of Eqs. (1), possibly supplemented with Eq.(2), 
is solved by the \ 2 method for unknown parameters e and 
(ef)- 

Care must be taken in calculating complex observa- 
tional amplitudes. If the time dependence of the observed 
quantity, O, is given in the form Ao sin(uit + <p) then the 
corresponding complex amplitude is (Ao sin ip, Ao cos ip). 
Comparing the empirical and theoretical / values, we have 
to make sure that in the nonadiabatic pulsation code the 
time dependence is e~ luJt . 

Effects of rotation are totally ignored in equations (1) 
and (2). The approximation is essential because only then 
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4,38 4,36 


4,34 


4,32 4,30 4,28 



Table 1. Parameters of selected models of S Ccti. 



'og^ff 

Fig. 1. The observational error box for S Ceti in the HR 
diagram. The evolutionary tracks were calculated for two 
values of the metallicity parameter Z. The lines of con- 
stant fundamental (n — 1) and the first overtone (n = 2) 
period equal to 0.161137 d. are shown. These lines will be 
discussed later in this section. 

unknown values of i and m may be absorbed into e. As 
long as the angular velocity of rotation is much lower 
than uj these equations provide a good approximation. 
Effects of rotation enter the expression for A x through 
terms cx f2 2 and the expression for M.^ through terms 
oc f2. Therefore use of Eq.(2) is less safe but, as we will 
see later, it is essential for mode identification in B-type 
stars. There is a need for improvement. Fortunately the 
two stars considered in this paper are slow rotators. 

In our applications we take uj from observations, the 
stellar radius, R, from the measured mean luminosity, 
L, and effective temperature, T c g, while the mass, M, is 
taken from evolutionary tracks. 

3. The monoperiodic j3 Cep star 5 Ceti 

3.1. Data 

We derive mean parameters of the star relying on 
Stromgren photometry from the catalogue of Stankov & 
Handler (2005) and on the Hipparcos parallax taking into 
account the Lutz-Kelker correction (Lutz & Kelker 1973). 
In Fig. 1 we show position of the star with its error box 
in the HR diagram. Table 1 presents 5 models which were 
chosen for analysis; all of them are consistent with the 
observations. The models correspond to the central value 
(Dl) and to four values close to the edges of the error box 
(D2-D5), and were fixed at Z = 0.02. 

We use here the amplitudes and phases in four 
Stromgren passbands from Kubiak & Seggewiss (1990). 
These photometric observations were done in 1986. We 
derived the corresponding radial velocity data from spec- 
tra taken in 1987 (Aerts et al. 1992). The data used in our 
analysis are given in Table 2. 





M/Mq 


log T cft 


logL 


log g 


Dl 


9.70 


4.3340 


3.930 


3.783 


D2 


9.00 


4.3229 


3.806 


3.830 


D3 


9.20 


4.3459 


3.790 


3.947 


D4 


10.20 


4.3232 


4.041 


3.650 


D5 


10.40 


4.3472 


4.032 


3.763 



Table 2. Amplitudes and phases in radial velocity and in 
four Stromgren passbands for monoperiodic (3 Cep star 6 
Cet with pulsation frequency 6.205875 c/d. 







A [km/s],[mag 


;] V? [rad] 








7.4323(723) 


1.6416(78) 






U 


0.0263(9) 


2.7758(156) 






V 


0.0143(9) 


2.7096(156) 






b 


0.0137(9) 


2.7291(156) 






y 


0.0128(9) 


2.7057(156) 




Table 3. Values of e, f and % 2 


for models of 5 Ceti 








El 


/r 


h 


x 1 






v = 6.20587 c 


/d, £ = 






Dl 


0.00520(5) 


0.00036(4) 


-9.70(18) 


2.56(12) 


1.03 


D2 


0.00570(5) 


0.00040(4) 


-9.13(15) 


2.36(10) 


0.74 


D3 


0.00645(6) 


0.00045(5) 


-8.37(14) 


2.07(10) 


0.90 


D4 


0.00435(5) 


0.00031(4) 


-1.12(27) 


3.05(17) 


1.55 


D5 


0.00491(5) 


0.00035(4) 


-1.01(22) 


2.69(14) 


1.28 



3.2. The I degree of the excited mode 

The only detected mode in this star has been already iden- 
tified as I — from spectroscopy by Aerts et al. (1992) 
and from photometry by Cugier et al. (1994). This identi- 
fication was confirmed by Cugier & Nowak (1997), who in- 
cluded observations in the UV range from IUE. However, 
in the photometric identification, the theoretical / val- 
ues were adopted. Naturally, extracting both £ and / 
from observations sets higher requirements on the data. 
Photometric data are insufficient. 

It is clear from comparison of the x 2 (^) dependence, 
shown in the two panels of Fig. 2, that radial velocity data 
are essential for determination of the spherical harmonic 
degree, I. We see in the upper plot that without radial 
velocity data, the discrimination between various even de- 
grees is not possible. The lower plot shows that the mode 
degree must be £ = 0. In addition, the lower luminosity 
models lead to a lower x 2 ■ In Table 3 we provide our solu- 
tions for the £ = identification. With this identification 
e is equal to the SR/R amplitude. 

3.3. Discrimination between two possible radial modes 

Having determined the degree of the excited mode we still 
have to consider its possible radial orders. In Fig. 1 we 
show the theoretical HR diagram with lines corresponding 
to a period of 0.16114 d for two different radial mode iden- 
tifications. Both fundamental and first overtone modes are 
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Fig. 2. The x 2 per degree of freedom as a function of I ob- 
tained from the fit of photometric amplitudes and phases 
without radial velocity data (upper panel) and with radial 
velocity data (lower panel). The five lines correspond to 
models from Table 1 which are consistent with the obser- 
vations. 

acceptable. Now, for each value of Z and n we have a one 
dimensional family of models to consider. 

We first check whether the x 2 behaviour may help us 
to discriminate between the two possible choices of n. The 
result shown in Fig. 3 is rather negative. The difference 
between the two choices is small, especially in the low 
temperature range, where we reach the lowest x 2 regard- 
less of the n values. On the same plot we show that the 
metallicity, [m/H], microturbulence, vt, and opacity in the 
atmosphere models do not matter very much for the x 2 
values. The metallicity affects the analysis in two ways. 
The first is through the evolutionary models, hence the 
surface gravity, the second is through the [m/H] parame- 
ter in the atmosphere. 

The lesson from the x 2 analysis is that the mode is def- 
initely the radial one and that there is some preference for 
lower effective temperature and luminosity within the er- 
ror box. However, the discrimination between mctallicitics 
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Fig. 3. The x 2 P er degree of freedom as a function of 
log T e s for two possible identifications of the radial mode 
order. The vertical lines show the allowed T c g range, con- 
sistent with the error box shown in Fig. 1. In the upper 
panel we show the effect of use of the new opacity distri- 
bution function (ODF) in the atmosphere and enhance- 
ment in O, Ne, Mg, Si, S, Ar, Ca and Ti (A+ncwODF), 
following Kurucz (2004). In the lower panel we compare 
the effect of different metallicity parameters, [m/H], and 
microturbulence velocities, v t . 



and two radial mode orders was not possible. Much more 
may be learnt from comparing the inferred nonadiabatic 
parameter, /, with model calculations. 

3.4. Inference from f -values and stability consideration 

As explained in the introduction, the calculated / values 
must be very sensitive to the iron abundance in the driv- 
ing zone. The effect is well seen in Fig. 4. We also see the 
sensitivity to the radial mode identification. Comparison 
of the upper and lower panel strongly suggests that the 
excited mode is the first overtone. This is the same identi- 
fication as proposed by Cugier et al. (1994) and supported 
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Fig. 4. Comparison of the empirical and theoretical val- 
ues of / assuming Z and T e g within the error box. The 
upper and lower panels refer to the fundamental and first 
overtone identification of the 5 Cet mode, respectively. 
The uncertainties in metallicity and effective temperature 
are included both in the empirical and theoretical values 
of /. The empirical values are taken from Table 3. The 
theoretical values were calculated along lines of constant 
period (see Fig. 1) with the step in mass of 0.2M©. In the 
upper panel, the masses at the top are 8.6, 8.8, 9.0 M & 
for Z — 0.015, 0.017, 0.020, respectively, whereas in the 
lower panel the masses at the top are 9.2, 9.4, 9.6 Mq for 
Z = 0.015, 0.017, 0.020, respectively. 



by Cugier & Nowak (1997). However there is a problem 
with this identification, because in the allowed ranges of 
T e ff and Z, the first overtone is stable, as shown in Fig. 5. 
The adopted measure of instability is 77, which varies be- 
tween — 1 and +1. For unstable modes we have 77 > 0. 
Though the mode is not far from the instability domain, 
we regard the discrepancy as serious because in the same 
range of parameters the fundamental mode is definitely 
unstable. This is a typical property of (3 Cep star mod- 
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too 7 "* 

Fig. 5. The ranges of T e g and Z allowed by consistency 
of empirical and calculated values of / assuming the first 
overtone identification. The lines of constant instability 
parameter, 77, are shown. Negative values of 77 correspond 
to stability. 



els, in which the lower frequency radial and neighbouring 
nonradial modes are preferentially unstable. 

We may think about various sources of the discrepancy 
in empirical as well as in theoretical values of /. The effect 
influencing the latter ones may be related to the opac- 
ity data. Fortunately, we have two independent sources of 
these data. Our results so far were obtained with the tab- 
ular opacities from the OPAL project (Iglesias & Rogers 
1996). As the alternative we have latest tables based on 
the Opacity Project (Seaton & Badnell 2004, Badncll ct 
al. 2005, Seaton 2005), hereafter OP. As for the relative 
heavy element abundance, the OPAL data use the GN93 
(Grevesse k Noels 1993) and the OP data use S92 (Seaton 
et al. 1994) solar mix. The two mixes differ only slightly. 
In particular, the S92 mix has a Fe abundance higher by 
2.5%. Different opacities lead only to small differences in 
the evolutionary tracks. The OP tracks are fainter by some 
0.02 dex. Surprisingly, a large difference was found in the 
/ values. Now, there is no clear discrimination between the 
fundamental and first overtone based on the comparison of 
the theoretical and empirical / values. Plots in Fig. 6 show 
that both identifications are at most marginally within 
with the error box. Note, however, that while for the fun- 
damental mode identification all acceptable models lie well 
within the instability domain (77 > 0), the first overtone 
identification places the models in the stability domain. 
This is why we prefer the fundamental mode identifica- 
tion, which implies, for instance, that luminosity is lower 
by 0.33 dex and mass is lower by 2 Mq than in the case 
of the first overtone identification. Still, improvements are 
needed for a convincing discrimination between the two 
possibilities. 
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Fig. 6. The same as Fig. 4 but calculated with the new 
OP opacity tables. In addition, lines of constant values 
of the instability parameter, 77, are shown. In the upper 
panel, the masses at the top are 8.8, 8.8, 9.0 M for Z 
= 0.017, 0.020, 0.023, respectively, whereas in the lower 
panel the masses at the top are 9.2, 9.4, 9.6, 9.8 M@ for 
Z = 0.015, 0.017, 0.020, 0.023, respectively. The step in 
mass is 0.2M Q . 
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Fig. 7. The observational error box for v Eri on the HR 
diagram. The evolutionary tracks for two values of the 
metallicity parameter Z are shown. 

Table 4. Parameters of selected models of v Eridani. 





M/Mq 


logT cff 


logL 


log 9 


Nl 


9.90 


4.3460 


3.940 


3.830 


N2 


9.20 


4.3344 


3.823 


3.869 


N3 


9.40 


4.3589 


3.792 


4.007 


N4 


10.40 


4.3331 


4.059 


3.681 


N5 


10.60 


4.3584 


4.044 


3.805 



We adopted the mean parameters of v Eri after 
Pamyatnykh, Handler & Dziembowski (2004), where de- 
tails about derivation are given. In Fig. 7, we show the 
star position in the HR diagram with the error box. 
Evolutionary tracks plotted in the figure were calcu- 
lated for two indicated Z values, assuming no convective 
overshooting. Seismic models of Pamyatnykh, Handler & 
Dziembowski (2004) were obtained at Z = 0.015. Table 
4 presents 5 models which were chosen for analysis, all of 
them are consistent with the observations. As in the case 
of 8 Ceti, these five models correspond to the central value 
(Nl) and to four values close to the edges of the error box 
(N2-N5), and were fixed at Z = 0.02. 



4. The multiperiodic (3 Cep star v Eridani 

This star has been the object of the recent photometric 
(Handler et al. 2004) and spectroscopic campaigns (Aerts 
et al. 2004, de Ridder et al. 2004). The number of detected 
modes is much larger than in any other star of this type. 
For several of them we have unambiguous mode identifica- 
tion, which was supported by seismic model constructions 
(Pamyatnykh, Handler & Dziembowski 2004, Ausseloos et 
al. 2004). 

Extracting / values from data for a multiperiodic ob- 
ject is clearly advantageous because we get independent 
constraints on mean stellar parameters, such as T e g or Z, 
from each mode. 



4.1. Data 

Photometric amplitudes and phases used in our anal- 
ysis of this star were kindly provided to us by M. 
Jerzykiewicz, who reanalyzed data from the multisite ob- 
servations carried out between October 2002 and February 
2003 (Handler et al. 2004). To determine the amplitudes 
and phases of the radial velocity, we calculated the first 
moment of the line profile of Silll 4553A using spectra 
obtained during the last spectroscopic campaign (Aerts 
et al. 2004). We stress that it is important to use con- 
temporaneous data in this case because amplitudes and 
phases do change at the level exceeding errors. It is also 
important that the same frequencies are used in analy- 
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Table 5. Amplitudes and phases in radial velocity and in 
four Stromgren passbands for four dominant frequencies 
in v Eri. 





A [km/s],[mag; 


Lp [rad] 




v x = 5.76326 


c/d 




22.017(116) 


5.473(5) 


It 


0.07345(20) 


0.616(3) 


V 


0.04103(14) 


0.590(3) 


y 


0.03685(13) 


0.579(4) 




v 2 = 5.65391 


c/d 


Vrad 


9.049(115) 


0.419(13) 


U 


0.03793(21) 


1.982(5) 


V 


0.02645(14) 


1.982(6) 


y 


0.02506(13) 


1.975(5) 




v-i = 5.62009 


c/d 


Kad 


8.272(112) 


0.129(14) 


U 


0.03464(21) 


1.681(6) 


V 


0.02388(15) 


1.691(6) 


y 


0.02267(14) 


1.684(6) 




v 4 = 5.63715 


c/d 


Kad 


7.703(105) 


3.033(14) 


u 


0.03217(22) 


4.537(7) 


V 


0.02244(15) 


4.537(7) 


y 


0.02105(14) 


4.544(7) 



sis of spectroscopic data. The latter was the reason why 
we could not rely on results of De Ridder et al. (2004), 
who analyzed data from the same campaigns. At the mo- 
ment, we have sufficiently accurate data only for the four 
dominant peaks. We have amplitudes and phases in three 
Stromgren passbands, uvy, and the radial velocity data. 
The data used in our analysis are listed in Table 5. 

4.2. The £ and f values 

Like in the case of 8 Cet, only after combining photometric 
and spectroscopic data, could we get a unique £ identifi- 
cation for the dominant modes in v Eri. Our £ degrees are 
the same as those of our predecessors. However, unlike 
them, we did not assume the / values. 

In Fig. 8, we compare % 2 values obtained with and 
without the radial velocity data for the dominant mode, 
v\, which has been identified as £ — (see e.g. De Ridder 
et al. 2004) on the base of photometric data alone but 
upon assuming theoretical /. We see in the upper part 
of the figure that our method without spectroscopic data 
does not allow for £ determination. In the case of combined 
data, the minimum of \ 2 is indeed very deep, however the 
value of x 2 is large (« 50). Also one gets a very large % 2 
(« 40) adopting / from theory and relying only on photo- 
metric data (De Ridder et al., 2004). For the three modes, 
v ii ^3, ^4, forming the close triplet, we confirm in Fig. 9 
the I = 1 identification. In Table 6 we provide results of 
our analysis for all four modes but only £ leading to the % 2 
minima. The minima for the triplet are significantly lower 
than for the £ = mode, but still larger than 1. All these 
modes have close frequencies thus they should have similar 
values of /. At frequencies around the fundamental radial 




1 2 3 „ 4 5 6 




I I i I i I . I . I i I i L 

1 2 3 „ 4 5 6 



Fig. 8. The \ 2 per degree of freedom as a function of £ ob- 
tained from the fit of photometric amplitudes and phases 
without radial velocity data (upper panel) and with radial 
velocity data (lower panel), for the dominant frequency, 
v\ . The five lines correspond to models from Table 1 which 
are consistent with the observations. 

mode and above, these values are nearly ^-independent, 
as long as £ is not large. It is so because for such modes 
pulsation is vertical in the layers contributing to /. Fig. 10 
shows that indeed all four / values determined from the 
data are close, which is reassuring. However, the large % 2 
are of concern. There are two possible explanations: either 
the model is inadequate or the errors are underestimated. 

Let us begin with the first possibility. Neglect of the 
effect of rotation in the model is fully justified because 
the star is a very slowly rotating object. The ratio of 
rotation to pulsation frequency is ~ 10~ 3 . More suspect 
may seem our use of linear relation between the displace- 
ment and the the relative flux perturbation which is only 
an approximation, but it seems justified. Consider the v\ 
mode, which was identified as radial and for which we 
determined the highest % 2 (still by far the lowest than 
with other identification). With such an identification we 
have |e| = |e| = 0.012 and |e/| = 0.11. The maximum 
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Fig. 9. The same as in the lower plot of Fig. 8 but for the close frequency triplet, 1/3, 1/4, V2 (£ = 1, m = —1, 0, +1) 



Table 6. Values of e, f and x 2 for models of v Eridani 







£l 


/r 


fx 


x J 


vi = 5.76326 c/d, I = 


Nl 


-0.0125(6) 


-0.0118(6) 


-9.2(4) 


0.8(4) 


48.4 


N2 


-0.0133(8) 


-0.0126(8) 


-8.8(5) 


0.7(5) 


79.6 


N3 


-0.0152(11) 


-0.0143(11) 


-8.1(5) 


0.6(5) 


112.0 


N4 


-0.0104(8) 


-0.0099(8) 


-10.6(7) 


0.9(7) 


128.4 


N5 


-0.0119(7) 


-0.0113(7) 


-9.6(4) 


0.8(4) 


61.8 




v-a = 


= 5.62009 c/d, I 


= 1, m = 


-1 




Nl 


0.0010(3) 


-0.0080(3) 


-9.0(3) 


0.1(3) 


6.81 


N2 


0.0011(3) 


-0.0085(3) 


-8.4(3) 


0.1(3) 


7.27 


N3 


0.0012(3) 


-0.0094(3) 


-7.5(2) 


0.1(2) 


4.15 


N4 


0.0009(3) 


-0.0067(3) 


-10.9(5) 


0.1(5) 


9.84 


N5 


0.0010(3) 


-0.0075(3) 


-9.6(4) 


0.1(4) 


8.22 




u A = 


= 5.63715 c/d, I 


— 1, m = 







Nl 


0.0008(2) 


0.0074(2) 


-9.0(2) 


0.5(2) 


3.02 


N2 


0.0009(2) 


0.0079(2) 


-8.4(2) 


0.5(2) 


3.40 


N3 


0.0010(1) 


0.0088(1) 


-7.5(1) 


0.5(1) 


1.46 


N4 


0.0007(2) 


0.0063(2) 


-10.9(3) 


0.6(3) 


4.95 


N5 


0.0008(2) 


0.0070(2) 


-9.6(3) 


0.6(3) 


3.93 




l>2 = 


= 5.65391 c/d, I 


= 1, m = 


1 




Nl 


0.0036(3) 


-0.0080(3) 


-9.1(3) 


0.1(3) 


7.07 


N2 


0.0038(3) 


-0.0086(3) 


-8.5(3) 


0.1(3) 


8.34 


N3 


0.0042(3) 


-0.0094(3) 


-7.5(2) 


0.1(2) 


4.06 


N4 


0.0030(3) 


-0.0068(3) 


-11.0(5) 


0.1(5) 


10.07 


N5 


0.003493) 


-0.0076(3) 


-9.7(3) 


0.1(3) 


7.49 



relative perturbation temperature is about 0.03, which 
is small. The most natural interpretation of the higher 
order peaks (harmonics and combinations) found e.g. by 
Handler et al. (2005) is a nonlinearity of the A x (e) depen- 
dence. The lowest-order nonlinear departure from Eq.(l) 
is cubic. There is a number of third peaks but there is 
none at The peak at 2i/i is seen with the amplitude 
1/16 of the v\ peak. We use these facts as an argument 
for the applicability of Eq.(l) in our case. 

The next questionable part of our method concerns us- 
age of equilibrium atmosphere models for evaluation of the 



flux derivatives in Eq.(l). Is the assumption that 8 J- and 
8r are constant within the atmosphere justified ? Are the 
tabular data sufficiently accurate ? After the discussion 
of Cugier et al. (1994), we are confident about a positive 
answer to the first question. We are less confident in the 
case of the second question but it seems unlikely that this 
is the cause of large x 2 because in the case of 8 Ceti the x 2 
minima were about 1. We used atmospheric models from 
the same source and the stellar parameters of these two 
stars are very close. 

Having rejected the inadequacy of the model, we turn 
to the error estimate. Let us first note that the main con- 
tribution to x 2 (above 90%) comes from photometric data 
and that, in the case of the v\ peak, we require a 7 times er- 
ror increase of the photometric measurement to get x ~ 1- 
The errors given in Table 6 could be underestimated. The 
differences in amplitudes derived from various data sets 
exceeds the errors by a factor 5. The problem was briefly 
discussed by Jerzykiewicz et al. (2005) who suggest only 
a factor of 2 error enhancement. A more sophisticated ap- 
proach to the error analysis is needed but it is beyond the 
scope of the present work. We proceed further relying on 
determined / values, being encouraged by the their con- 
sistency and reasonably low errors as well by the clear i 
discrimination. 

In our comparison of empirical and theoretical values 
of /, we will consider only the I — mode because the 
errors are the lowest in that case. It is also safer than us- 
ing the average from the four modes because there may 
be a noticeable effect of difference in frequencies. Values 
of / for modes at more distant frequencies would be much 
more useful for estimating the j(y) dependence, which is 
of interest. Unfortunately, we did not succeed in a mean- 
ingful determination of / for any high frequency mode in 
v Eri. 

Comparison of the two panels of Fig. 11 shows clearly 
that the fundamental mode identification is preferred. The 
same identification was obtained by Pamyatnykh, Handler 
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Fig. 10. Empirical / values for v\, (I = 0), and the 
close triplet ^2,^3,^4, [t = 1) for logT e g = 4.355 and 
log L/L = 3.969. 

& Dziembowski (2004) and Ausseloos et al. (2004) who 
constructed seismic models v Eri based on a simultaneous 
fit of the radial and dipole mode frequencies measured in 
the star. 

As in the case of the 8 Cet model, we find that opac- 
ity significantly influences the theoretical / but, regard- 
less of the source of the opacity data, the fundamental ra- 
dial mode is the only acceptable identification. Comparing 
Fig. 12 with the upper panel of Fig. 11, we see somewhat 
better agreement between the empirical and theoretical 
values when the OP opacities are used. However, the in- 
ference on the star parameter is significantly different. For 
instance, using OP data we get the mass higher by about 
O.4M , Z lower by about 0.0015, and logT cff higher by 
0.01. The above conclusions are based on stellar models 
constrained only by mean parameters from photometry 
and the radial mode frequency, v\. We will see that our 
new seismic models lead us to somewhat different results. 

4.3. New seismic models 

We constructed seismic models using both OPAL and OP 
data. Like Pamyatnykh, Handler & Dziembowski (2004), 
the construction uses frequencies of the I = 0,pi, I = l,gi 
and £ — l,pi modes. Models calculated with OPAL are 
slightly different than those of Pamyatnykh, Handler & 
Dziembowski (2004) due to a higher accuracy demand on 
the frequency fit. The model parameters and the / values 
for i = 0,pi are listed in Table 7. For each of the opacities, 
we provide data for the models obtained with two values 
of the overshooting parameter, a ov — and 0.1. In the 
seismic models, the influence of a ov on / enters indirectly 
through T c g and Z. We may see that the parameters of the 
models are significantly affected by the choice of the opac- 
ity data. Models calculated with the OP data are cooler 
by about 0.015 dex. The model calculated with a ov = 0.1 



2.0 

1.5 

1.0 

0.5 |- 

0.0 
-0.5 
-1.0 h 
-1.5 



all unstable 




\ 



M, log7~ all 



Z=0.015 



Z =°° 17 Z=0.020 



-10.5 -10.0 -9.5 -9.0 -8.5 -8.0 -7.5 

L 



4.5 
4.0 
3.5 

f 3.0 

'1 

2.5 - 
2.0 - 
1.5 
1.0 
0.5 
0.0 I- 



all stable 




\M,iogr„, 

Z=0.015 -7 nn-iv 

Z=0.017 z = o.020 




-12.0 -11.5 -11.0 -10.5 -10.0 -9.5 
L 



-9.0 



Fig. 11. Theoretical and empirical / for the radial mode 
(v\) assuming various stellar parameters within the error 
box. The upper panel assumes fundamental and the lower 
one the first overtone mode identification. The empirical 
values are taken from Table 6. Theoretical / were obtained 
adopting the three indicated values of the metallicity pa- 
rameter, Z , and the OPAL opacities. In the upper panel, 
the masses at the top are 9.0, 9.0, 9.2 M for Z = 0.015, 
0.017, 0.020, respectively, whereas in the lower panel the 
masses at the top are 9.6, 9.8, 10.0 M for Z = 0.015, 
0.017, 0.020, respectively. The step in mass is 0.2M Q . 



is in fact somewhat outside the allowed range of T e g . The 
OP models have higher Z by about 0.003. We see that 
largest differences between the OPAL and OP models are 
in the imaginary part of /. Indeed the value of / is a very 
sensitive probe of the stellar opacities. 

How the / values from the four seismic models compare 
with observations is shown in Fig. 13. The correspond- 
ing empirical / are also given in Table 8. For none of 
the opacity data is the agreement fully satisfactory. The 
truth seems to be somewhere in between, if a moderate 
overshooting is allowed. Relying only on OPAL opacities, 
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Fig. 12. The same as in the upper panel of Fig. 11 but 
with the theoretical / values calculated with the OP opac- 
ities. Models lying to the right of the line marked with 
r\ = 0.0 are unstable. The masses at the top are 8.8, 9.0, 
9.2 M Q for Z = 0.015, 0.017, 0.020, respectively. The step 
in mass is 0.2M Q . 

Table 7. Seismic models of v Eri calculated with the 
OPAL and OP opacities. The / values are given for 
t = 0,Pi- 





£R 


£l 


fR 


ft 


x 2 




v\ 


= 5.76326 


c/d, 1 = 






SI 


-0.0126(6) 


-0.0120(6) 


-9.2(3) 


0.8(3) 


40.2 


S2 


-0.0128(7) 


-0.0121(7) 


-9.1(4) 


0.8(4) 


55.8 


S3 


-0.0126(7) 


-0.0120(7) 


-9.2(4) 


0.8(4) 


56.4 


S4 


-0.0127(7) 


-0.0121(7) 


-9.2(4) 


0.8(4) 


61.9 



we may reach the agreement by increasing a ov to about 
0.2. However, a further increase would leave the instability 
range. 

5. Conclusions 

Using amplitudes and phases from simultaneous photo- 
metric and spectroscopic observations of two j3 Cep stars, 
we derived the spherical harmonic degree, £, of excited 
modes and the complex parameter, /, which may be di- 
rectly compared with the nonadiabatic theory of stellar 
oscillations. The / parameter determines bolometric flux 






M/Mq 


CXov 


z 


logT cff 


logL 


fR 


// 


OPAL 


SI 


9.808 


0.0 


0.0155 


4.3530 


3.959 


-8.56 


-0.63 


S2 


9.230 


0.1 


0.0145 


4.3424 


3.903 


-8.82 


-0.07 


OP 


S3 


9.590 


0.0 


0.0185 


4.3396 


3.902 


-8.22 


1.20 


S4 


9.020 


0.1 


0.0175 


4.3284 


3.843 


-8.70 


1.70 



Table 8. Values of e, f and % 2 f° r seismic models of v 
Eridani for the dominant frequency, v\. 



Fig. 13. The f values for fundamental radial mode of the 
seismic models constructed using opacity data from two 
sources (OPAL and OP) and using two values of the over- 
shooting parameter (a ov = 0.0 and 0.1). The empirical 
f values were obtained assuming atmospheric parameters 
corresponding to the seismic models. 



amplitude at specified surface displacement and it consti- 
tutes a new seismic probe of stellar outer layers. 

One of the stars is a monoperiodic pulsator 5 Cet, the 
other is a multipcriodic pulsator v Eri. In the latter star, 
we analyzed data on four dominant peaks. We found that 
unique determination of mode degree requires both pho- 
tometric and spectroscopic data. In our analysis we relied 
on the use of static atmospheric models of Kurucz (2004). 

Empirical values of / were compared with the theo- 
retical values for relevant stellar models. The models em- 
ployed the OPAL (Iglesias & Rogers 1996) and the OP 
tabular data (Seaton & Badnell 2004, Badnell et al. 2005, 
Scaton 2005). The evolutionary tracks are not much af- 
fected by the choice of the data but significant differences 
are seen in the resulting /. 

We confirmed the I — identification for the only 
mode detected in 5 Cet but, unlike our predecessors, we 
think that the mode is more likely fundamental than the 
first overtone. The reason is that the latter mode is stable 
over the allowed range of mean star parameters. If indeed 
the mode is fundamental then approximate consistency of 
the / values is possible only when the OP opacities are 
used. 

Our identification of the four dominant peaks in 
the v Eri oscillation spectrum is the same as that 
of Pamyatnykh, Handler & Dziembowski (2004) and 
Ausseloos et al. (2004). Consistency between the theo- 
retical and empirical values of / was achieved both with 
OPAL and OP data, but at significantly different stellar 
parameters. 

For none of the two opacity data did we achieve full 
agreement between empirical / and those derived from 
seismic models fitting three dominant centroid frequencies 
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in the oscillation spectrum. As one can see, a perfect fit 
would be possible at some intermediate opacities assuming 
a small overshooting parameter, a ov . 

There are things to improve in our models of Cep 
stars. As Pamyatnykh, Handler & Dziembowski (2004) 
pointed out, the standard models do not account for ex- 
citation of modes at very low and very high frequen- 
cies in v Eri. We also cannot be fully satisfied with our 
fit of empirical and theoretical /. These two problems 
may be related. The solution to the first problem pro- 
posed by Pamyatnykh, Handler & Dziembowski (2004) 
was a significant enhancement of the iron abundance in 
the driving zone caused by the selective radiation pres- 
sure. Undoubtedly, this will be reflected in the / values. 
We plan to investigate the effect when we have physically 
justified models of the enhancement. 

Ideal seismic stellar models should account, within the 
errors, not only for all measured frequencies but also for 
the / values. Such models will provide an ultimate test of 
our knowledge of the stellar structure from the atmosphere 
to the deep interior. We are not at this point yet. We 
believe that the most important result of the present work 
is that our new seismic tool provides a stringent probe of 
stellar opacities. 

Acknowledgements. We are grateful to Mike Jerzykiewicz for 
kindly providing data on photometric amplitudes and phases 
as well as for illuminating discussion about possible sources of 
errors. We thank Mike Seaton for detailed discussion of how to 
use his new opacity data, and Conny Aerts for making spectra 
of S Cet and v Eri accessible. We also thank the anonymous 
referee for stimulating comments. JDD thanks the Foundation 
for Polish Science for supporting her stay in the Copernicus 
Astronomical Center. The work was partially supported by 
Polish MNil grant No. 1 P03D 021 28. 



Handler, G., Shobbrook, R. R., Jerzykiewicz, M. et al., 2004, 

MNRAS, 347, 454 
Iglesias, C. A., Rogers, F. J., 1996, ApJ 464, 943 
Jerzykiewicz, M., Handler, G., Shobbrook, R. R., et al., 2005, 

MNRAS, 360, 619 
Kolaczkowski, Z., Pigulski, A., Soszyhski, I., et al., 2004, in 

Variable Stars in the Local Group, Proc. IAU Coll. 193, 

eds. D. W. Kurtz and K. R. Pollard, ASP Conf. Ser., Vol. 

310, 225 

Kubiak, M., Seggewiss, W., 1990, Acta Astron. 40, 85 
Kurucz, R. L., 2004, http://kurucz.harvard.edu 
Lutz, T. E., Kelker, D. H., 1973, PASP 85, 573 
Pamyatnykh, A. A., Handler, G., Dziembowski, W. A., 2004, 

MNRAS 350, 1022 
de Ridder, J., Telting, J. H., Balona, L. A., et al., 2004, 

MNRAS 351, 324 

Seaton, M. J., 2005, MNRAS, in press, (astro-ph/04110101 
Seaton, M. J., Badnell, N. R., 2004, MNRAS 354, 457 
Seaton, M. J., Yu Yan, Mihalas, D., Pradhan, A. K., 1994, 

MNRAS, 266, 805 
Stankov, A., Handler G., 2005, ApJS, 158, 193 



References 

Ausseloos, M., Scuflaire, R., Thoul, A., Aerts, C, 2004, 

MNRAS 355, 352 
Aerts, C, de Cat, P., Handler, C, et al., 2004, MNRAS 347, 

463 

Aerts, C, de Pauw, M., Waelkens, C, 1992, A&A 266, 294 
Badnell, N. R., Baut ista, M. A., Butl er, K., et al., 2005, 

MNRAS, in press (astro-ph/0410744 ) 
Claret, A., 2000, A&A 363, 1081 

Cugier, H., Dziembowski, W. A., Pamyatnykh, A. A., 1994, 

A&A 291, 143 
Cugier, H., Nowak, D., 1997, A&A 326, 620 
Daszyhska-Daszkiewicz, J., Dziembowski, W. A., Pamyatnykh, 

A. A., Goupil, M-J. 2002, A&A 392, 151 
Daszyhska-Daszkiewicz J., Dziembowski W. A., Pamyatnykh 

A. A., 2003, A&A 407, 999 
Daszyhska-Daszkiewicz, J., Dziembowski, W. A., Pamyatnykh, 

A. A., Breger, M., Zima, W., 2004, Proc. IAU Symp. 

224, eds. J. Zverko, W. W. Weiss, J. Ziznovsky and S. J. 

Adelman, p. 853 
Grevesse, N., Noels, A., 1993, in Origin and Evolution of the 

Elements, eds. N. Pratzo, E. Vangioni-Flam and M. Casse, 

Cambridge Univ. Press, p. 15 



